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HsTXl toxin, from the scorpion Heterometrus spinnifer, is a 34-residue, C-terminally amidated peptide 
cross-linked by four disulfide bridges. Here we describe new HsTXl analogues with an Ala, Phe, Val or Abu 
substitution at position 14. Complexes of HsTXl with the voltage-gated potassium channels Kvl. 3 and 
Kvl.l were created using docking and molecular dynamics simulations, then umbrella sampling 
simulations were performed to construct the potential of mean force (PMF) of the ligand and calculate the 
corresponding binding free energy for the most stable configuration. The PMF method predicted that the 
R14A mutation in HsTXl would yield a > 2 kcal/mol gain for the Kvl.3/Kvl.l selectivity free energy 
relative to the wild-type peptide. Functional assays confirmed the predicted selectivity gain for 
HsTXl [R14A] and HsTXl [R14Abu], with an affinity for Kvl. 3 in the low picomolar range and a selectivity 
of more than 2,000-fold for Kvl. 3 over Kvl.l. This remarkable potency and selectivity for Kvl. 3, which is 
significantly up-regulated in activated effector memory cells in humans, suggest that these analogues 
represent valuable leads in the development of therapeutics for autoimmune diseases. 



The voltage-gated potassium channel Kvl. 3 is significantly up-regulated in activated effector memory (Tem) 
cells in humans^'^. As a consequence, Kvl. 3 blockers constitute valuable new therapeutic leads for the 
treatment of autoimmune diseases mediated by Tem cells, such as multiple sclerosis and rheumatoid 
arthritis^"^. 

HsTXl toxin is a 34-residue, C-terminally amidated peptide from the scorpion Heterometrus spinnifer, which 
is cross-linked by four disulfide bridges^. Its amino acid sequence and the locations of the four disulfide bridges 
are shown in Fig. 1. The solution structures of the synthetic toxin^ and a chimera consisting of the N-terminal half 
of the closely-related scorpion toxin maurotoxin and the C-terminal half of HsTXl ^ were found to be very similar 
to the canonical fold adopted by related scorpion toxins that contain only three disulfides, which consists of an N- 
terminal helix structure connected to a C-terminal two-stranded antiparallel p-sheet. 

HsTXl is a potent blocker of potassium channels. Lebrun et al.^ found that the peptide inhibited rat Kvl. 3 
channels with an IC50 of ca 12 pM and that it did not compete with ^^^I-apamin for binding to rat brain 
synaptosomal membranes, although it did compete efficiently with ^^^I-kaliotoxin for binding to voltage-gated 
K^ channels on the same preparation (IC50 ca 1 pM). It is thus a more potent Kvl. 3 channel blocker than 
maurotoxin and far more specific^. Subsequent studies confirmed its marked preference for Kvl. 3 over Kvl.l, 
Kvl. 2 and Kca3.1 channels^. This selectivity for Kvl. 3 makes HsTXl a potentially attractive candidate for the 
treatment of autoimmune diseases such as multiple sclerosis and rheumatoid arthritis, as it has been shown that 
blockade of this channel in self-reactive Tem cells is effective in preventing the tissue damage associated with these 
conditions'"^. Indeed, a considerable effort has been devoted to developing Kvl. 3 -selective analogues of the sea 
anemone peptide ShK as therapeutics for multiple sclerosis^'^"'" and one of these analogues has recently entered 
phase 1 clinical trials '\ 
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Figure 1 | Structure of HsTXl. (A) Amino acid sequence of HsTXl 
(UniProtKB/Swiss-Prot P59867; KAX63_HETSP) showing the positions 
of the four disulfide bridges^ (B) Simulated structure of HsTXl [R14A] 
(green) superposed on the solution structure of HsTXl ^ (transparent grey) 
to show that the R14A mutation does not cause any structural changes. 
Both structures are obtained from equilibrated MD simulations. The side 
chains of residue 14 and Lys23 in both structures are shown (HsTXl in 
grey, HsTXl [R14A] in colour). 

In an effort to understand the molecular basis for the potency and 
selectivity of HsTXl for Kvl.3, we recently undertook a computa- 
tional study of its interaction with Kvl.l, 1.2 and 1.3^^. Accurate 
models of Kvl.x-HsTXl complexes were created using docking 
and molecular dynamics simulations. For each complex, the binding 
free energy of HsTXl was determined from the potential of mean 
force calculations, with good agreement being obtained between the 
computed and experimental binding free energies. Comparison of 
the binding modes of HsTXl with Kvl.l and Kvl.3 revealed that the 
lower affinity of HsTXl for Kvl.l is due to its inability to come close 
to the pore domain, which prevents the pore-inserting lysine from 
making proper contacts with the tyrosine carbonyls in the selectivity 
filter of the channel. 

In this study we have utilized these models to design analogues of 
HsTXl with mutations at position 14, which resulted in even greater 
selectivity for Kvl.3 over Kvl.l and other channels. Because the R14 
mutations perturb the binding mode of HsTXl, the change in bind- 
ing free energy associated with these mutations may not be reliably 
calculated using path-independent methods such as free energy per- 
turbation. We therefore performed umbrella sampling molecular 
dynamics (MD) simulations and determined the binding free ener- 
gies of the peptide and HsTXl [R14A] from potential of mean force 
(PMF) calculations. The PMF method predicts that the R14A muta- 
tion in HsTXl will yield a > 2 kcal/mol gain in the Kvl.3/Kvl.l 
selectivity free energy relative to the wild-type peptide. Functional 
assays confirm the predicted selectivity gain for HsTXl [R14A] and 
other R14 analogues, and suggest that they will be valuable leads in 
the development of therapeutics for autoimmune diseases. 

Results 

Interactions of HsTXl [R14A] with Kvl.3 and Kvl.l. The desire to 
modify HsTXl at a site distant from its key pharmacophore led us to 
utilize docking simulations with both Kvl.l and Kvl.3 to predict 
residues that would improve selectivity for Kvl.3 over Kvl.l. 
Comparing the binding modes of HsTXl with Kvl.l and Kvl.3^^, 
we observed that the R14 side chain made an ionic interaction with 
the E353 side chain in the turret region of Kvl.l but had no contacts 



with Kvl.3 residues. This suggested that replacing R14 with a neutral 
residue could reduce the toxin's affinity for Kvl.l without affecting 
its potency for Kvl .3. We therefore explored in this study the effect of 
replacing R14 in HsTXl with alternative side chains. Snapshots of the 
Kvl.x-HsTXl[R14A] complexes are shown in Fig. 2. To facilitate 
comparison of the binding mode of this analogue with that of the 
wild-type toxin^^, we have superimposed the Kvl.x-HsTXl structure 
on that of HsTXl [R14A] after aligning the two structures. It is 
evident that, although there are no large rotations, the structure of 
HsTXl [R14A] is shghtly shifted relative to HsTXl. The effects of 
these changes on the binding modes are quantified in Table 1, which 
lists the average pair distances for the strongly interacting residues. In 
Kvl.l, the lost R14 contact is replaced by the N26-Y379 H-bond. In 
addition, R33 switches from E353 to E351, making a weaker contact. 
In Kvl.3, there are still five contacts but two of them (Y21 and N26) 
are replaced by R27 and K30. From the comparison of the binding 
modes, we expect the R14A mutation to reduce the affinity of the 
toxin to Kvl.l but not to Kvl.3, thus improving the Kvl.3/Kvl.l 
selectivity margin. To calculate by how much, we constructed the 
PMFs for the Kvl.x-HsTXl [R14A] complexes. 

Umbrella sampling MD simulations were performed as described 
in our earlier work on binding of toxin peptides to potassium 
channels^'^'^^'^^. Inspection of the backbone RMSD values of 
HsTXl [R14A] in each umbrella window showed that they exhibit 
a similar behaviour to those of HsTXl Namely, binding suppresses 
the RMSD values relative to the bulk, but they return to the bulk 
values after the toxin unbinds (data not shown). Thus, no artificial 
deformation of the peptide occurs during the PMF calculations. 
Convergence of the PMFs is demonstrated in supplementary Fig. 
SI. For Kvl.3, the PMF profile of HsTXl [R14A] and its convergence 
characteristics are very similar to those of HsTXl (Figure 5 in Rashid 
& Kuyucak^'). Despite weaker binding of HsTXl [R14A] to Kvl.l 
relative to HsTXl, its PMF converges much faster. To understand 
this feature, it is necessary to examine how the pair distances change 
as the toxin is pulled out. These are shown in supplementary Fig. S2 
for both complexes. Comparison of the pair distances for Kvl.3- 
HsTXl[R14A] with those of Kvl.3-HsTXl (supplementary Fig. 
S2) shows that the pair distances exhibit a very similar behaviour 
in both complexes. In contrast, there is a notable difference between 
the pair distances in the Kvl.l complexes. The large fluctuations in 
the R4 and R33 pair distances seen in the unbinding of HsTXl occur 
only for the R4 pair distance in HsTXl [R14A] , while R33 dissociates 
in the first step of pulling and never makes a contact again. The 
reduction in the side chain fluctuations is presumably responsible 
for the faster convergence of the Kvl.l -HsTXl [R14A]. More impor- 
tantly, the R14A mutation caused further destabilization of the 
binding mode, leading to a rapid detachment of the R33 contact. 
While this is compensated to a degree by the formation of the 
N26-Y379 H-bond, we expect the overall effect of the R14A muta- 
tion on the binding free energy of the toxin on Kvl.l to be more than 
the 2 kcal/mol value associated with the typical loss of a charge 
contact. 

The PMFs for HsTXl [R14A] are compared to those of HsTXl in 
Fig. 3, and the absolute binding free energies obtained from the 
integration of the PMFs are listed in Table 2. The values of R used 
in the calculation of the binding constants from eq. 1 are 0.72 and 
0.70 A, respectively, for Kvl.l and Kvl.3. The PMFs for Kvl.3 have 
almost the same well depth (Fig. 3) and the corresponding binding 
free energies are also very similar (Table 2). In contrast, the R14A 
mutation reduced the well depth in Kvl.l by 3 kcal/mol and 
increased the binding free energy of the toxin by 2.7 kcal/mol. 
Both of these results are in line with the qualitative expectations 
gathered from comparison of the binding modes, namely, that the 
R14A mutation should not affect the toxin's affinity to Kvl.3 much 
but should reduce its affinity for Kvl.l. In summary, the R14A 
mutation has increased the Kvl.3/Kvl.l selectivity free energy from 
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Figure 2 | Interaction of HsTXl [R14A] with potassium channels. Comparison of the binding modes of HsTXl [R14A] (green with coloured side chains) 
and HsTXl (transparent grey) in complex with the Kvl.l and Kvl.3 channels. The R14A mutation causes a small shift in the binding modes of the 
toxin in both complexes, resulting in some rearrangement of the contact pairs, as documented in the text and Table 1. 



3.9 kcal/mol in HsTXl to 6.2 kcal/mol in HsTXl [R14A]. This is 
sufficient to render HsTXl [R14A] effectively inactive on Kvl.l. 

Peptide synthesis and conformation. In order to confirm the 
prediction of enhanced selectivity for Kvl.3 over Kvl.l associated 
with substitution at Argl4, we replaced this residue in the native 
sequence with Ala, Abu, Phe or Val. The peptides were assembled 
using standard Fmoc-tBu solid-phase peptide synthesis. The crude 
products were oxidized using the glutathione-mediated oxidative 
folding conditions that have been used successfully for many other 
ShK analogues^'^^'^l As shown in Fig. 4, HsTXl [R14A] folded rapidly 
to a single major product, resulting in the typical pattern of a major 
earlier- eluting peak by RP-HPLC followed by later-eluting misfolded 
species and side-products. Each HsTXl analogue was purified to 
homogeneity by preparative RP-HPLC, as illustrated in Fig. 4 for 
HsTXl [R14A]. Each peptide had the correct mass by ESI-MS (data 
not shown), demonstrating that the four disulfide bonds had been 
formed. The yield was —40% of theory based upon the starting resin. 

CD spectra were acquired on both HsTXl and HsTXl [R14A] in 
dilute acetate buffer. Both peptides exhibited minima at around 
221 nm and maxima at around 195 nm (data not shown), very 



similar to those reported by Lebrun et al.^ and indicative of a folded 
peptide containing both a-helical and P-sheet secondary structure. 

One-dimensionaPH NMR spectra of HsTXl [R14A] showed sharp 
and well- dispersed resonances similar to those of wild-type HsTXl 
(Fig. 5), indicating that the R14A mutation did not cause any signifi- 
cant perturbation of the native structure. To further verify the struc- 
tural similarity of wild-type HsTXl and HsTXl [R14A], chemical 
shift differences from random coil values for backbone amide and 
H*^ resonances were plotted (supplementary Table SI and Fig. S3). 
This shows that the pattern of chemical shift deviations from random 
coil for HsTXl [R14A] closely resembles that of wild-type HsTXl, 
and confirms that the three-dimensional structures are basically the 
same. A comparison of the wild-type structure^ with the simulated 
structure of HsTXl [R14A] is shown in Fig. IB. 

Potenq^ and selectivity of potassium channel blockade. We used 
whole-cell patch-clamp to experimentally validate the modelling 
results. We patch-clamped mouse fibroblasts stably expressing 
either Kvl.3 or Kvl.l channels. Perfusion of HsTXl induced a 
block of both Kvl.3 and Kvl.l currents with IC50S of 29 ± 3 pM 
andll,330± 1,329 pM, respectively (Fig. 6 A, B). As predicted by the 
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Table 1 


Comparison of the strongly interacting pair distances 


in the HsTXl -Kvl .x compi 


exes with those 


in the HsTXl [R14A] -Kvl .x 


complexes 












HsTXl 




Kvl.l 


MD average 


HsTXl [R14A] 


Kvl.l 


MD average 


R4-N2 




E353-02(C) 


3.0 ± 0.5 


R4-N2 


E353-02(C) 


2.7 ±0.3 
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E353-Oi(B) 


2.9 ± 0.3 
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Y375-0(C) 
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K23-Ni 
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nces obtained from the MD simulations are listed in columns 3 and 6 (in units of A). Subscripts refer to the side 


chain atoms and the monomer identity is indicated in parentheses. 
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Figure 3 | PMFs for the unbinding of HsTXl [R14A] and HsTXl from 
Kvl.l and Kvl. 3. The R14A mutation has httle effect on the binding of 
HsTXl to Kvl. 3 but reduces its affinity for Kvl.l substantially. 



modelling data, HsTXl [R14A] retained its high affinity for Kvl. 3 
(IC50 = 45 ± 3 pM) and exhibited more than 2,000-fold selectivity 
for Kvl. 3 over Kvl.l with less than 20% Kvl.l current blocked at a 
dose of 100 nM (Fig. 6 A, B). Substitution of R14 with Abu generated 
the peptide HsTXl [R14Abu], which displayed a similar affinity and 
selectivity for Kvl. 3 over Kvl.l with IC50 ~ 50 pM on Kvl. 3 and > 
100 nM on Kvl.l (data not shown). 

Inhibitory effects on lymphocyte proHferation. Since CCR7" Tem 
and NA-NK lymphocytes are exquisitely sensitive to Kvl. 3 channel 
block^'^'^^, we tested the effects of HsTXl and the analogues 
HsTXl [R14A], HsTXl [R14Abu], HsTXl [R14F] and HsTXl [R14V] 
on the proliferation of rat CCR7" Tem cells, rat CCR7^ naive/TcM 
cells, and human A-NK and NA-NK cells. As a positive control, we 
used the well- characterized Kvl. 3 channel blocker ShK-186^'^. All 
peptides inhibited the proliferation of Ova-GFP Tem lymphocytes 
with IC50S ~ 1 nM (Fig. 6C). In contrast, a concentration of 
100 nM of these peptides was necessary to induce a —50% block in 
proliferation of splenic naive/TcM lymphocytes (Fig. 6C). ShK-186, 
HsTXl, and its analogues inhibited the mitogen-induced prolife- 
ration of NA-NK lymphocytes with IC50S of 5-15 nM but had no 
significant effect on the proliferation of A-NK cells (Fig. 6D). 

Peptide stability. Since HsTXl [R14A] and other R14 analogues 
displayed high potency and selectivity for Kvl. 3 and a corres- 
ponding ability to suppress lymphocyte proliferation, we also 
assessed their resistance to proteolysis in vitro. HsTXl [R14A] 
displayed no sign of degradation in presence of pepsin (Fig. 7A), 
presumably because of the lack of sequence-specific cleavage sites 
for this enzyme. In contrast limited digestion (—20%) was observed 



Table 2 | Comparison of binding free energies for the HsTXl - 
Kvl .X and HsTXl [Rl 4A] -Kvl .x complexes 



Complex 

Kvl.l -HsTXl 

Kvl. 1 -HsTXl [R14A] 

Kvl.3-HsTXl 

Kvl.3-HsTXl[R14A] 



AG,, 



3// 



Gb (PMF) 



-13.5 ± 0.3 
-10.5 ± 0.3 
-17.2 ± 0.3 
-17.1 ± 0.4 



-10.1 : 
-7.4 : 
-14.0 : 
-13.6 : 



0.3 
0.3 
0.3 
0.4 



-1 1.1 ± 0.1 
-14.9 ± 0.2 



The relative binding free energies obtained from the well depth in the PMFs are shown in the second 
column. The standard binding free energies determined from the PMFs using eq 2 (third column) 
are compared to the experimental values in the last column. Errors in the binding free energies are 
estimated from the block data analysis of the PMF data. All energies are in kcal/mol. 
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Figure 4 | Analysis of HsTXl and HsTXl[R14A]. RP-HPLC profiles for the crude refolded HsTXl (A) and HsTXl[R14A] (C), and purified 
refolded HsTXl (B) and HsTXl [R14A] (D). The gradient was from 10-40% B in 30 min at a flow rate of 1 mL/min. Buffer A was 0.05% TFA in H2O and 
Buffer B was 0.05% TFA in MeCN. The A220 was monitored for peak elution. 



HsTX1[R14A] 




''H chemical shift (ppm) 

Figure 5 | Conformation of HsTXl [R14A] . Comparison of ID^H NMR spectra of wild-type HsTXl (bottom) and HsTXl [R14A] , both in water at pH 
4.0, and acquired on a Bruker Avance 600 MHz spectrometer at 35°C. The amide-aromatic regions of both spectra are expanded 2.5-fold vertically 
relative to the aliphatic regions. 
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Figure 6 | Functional analyses of HsTXl[R14A]. (A) Whole-cell Kvl.3 (left) and Kvl.l (right) currents measured by patch -clamp in stably transfected 
L929 fibroblasts before (control) and after perfusion of HsTx-1 (top panels) or of HsTXl [R14A] (bottom panels). (B) Dose-response inhibition of Kvl.3 
(left) and Kvl.l (right) currents by HsTXl ( O and dashed line) and HsTXl [R14A] ( ■ ) fitted to a Hill equation (N = 3 cells per concentration). 
(C) Effects ofShK-186 (□ and grey line), HsTXl (O and dashed line), and HsTXl [R14A] (■) on the proliferation of rat Ova-GFP Tem lymphocytes 
(left) and of rat splenic T lymphocytes [mainly naive/TcM cells] (right) measured ex vivo by the incorporation of [^H] thymidine in the DNA of dividing 
cells (N = 3). (D) Effects ofShK-186 (□ and grey line), HsTXl (O and dashed line), and HsTXl [R14A] (■) on the proliferation of human NA-NK 
(left) and A-NK (right) lymphocytes (N = 3). 



by oc-chymotrypsin and trypsin (Fig. 7B,C) after 4 h at 37°C even 
though HsTXl [R14A] contains numerous basic and aliphatic 
residues that represent potential sites of cleavage. 

We also investigated whether these peptides caused any haemolysis 
in human erythrocytes. As shown in Fig. 7D, none of the peptides 
displayed haemolytic activity when tested at doses of 100 nM. 

Discussion 

We have created several analogues of the scorpion toxin HsTXl in 
which Argl4 is replaced with non-polar aliphatic or aromatic resi- 
dues. HsTXl [R14A] retains high affinity for Kvl.3 (IC50 = 45 ± 



3 pM) and exhibits more than 2,000-fold selectivity for Kvl.3 over 
Kvl.l. As selectivity for Kvl.3 over Kvl.l is an important criterion 
for potential therapeutic applications in the treatment of auto- 
immune disease in humans^'^, this peptide shows considerable prom- 
ise as a new potential therapeutic lead. 

With four disulfide bridges, HsTXl has the potential to generate 
105 distinct disulfide isomers during oxidative refolding in vitro. Our 
results for both the native peptide and HsTXl [R14A] show, however, 
that refolding is highly efficient and yields a single major isomer in 
good yield. In fact the presence of four disulfide bridges proves to be a 
significant benefit as the folded peptide is highly stable over a range of 
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Figure 7 | Stability of HsTXl[R14A]. Reversed-phase HPLC analysis after treatment with pepsin (A), a-chymotrypsin (B) and trypsin (C). 
(D) Haemolysis assay; comparison of the haemolytic ability of HsTXl and its analogues to haemolysis induced by Triton X-100 indicates that 
neither HsTXl nor its analogues induces significant haemolysis. Data shown as average ± SEM (N=3). 



pH and temperature in aqueous solution (our unpublished results). 
Consistent with this stability, HsTXl [R14A] is resistant to proteo- 
lysis by pepsin and showed only limited susceptibility to trypsin and 
chymotrypsin. There is an increasing acceptance of peptides as drug 
candidates'^ and an analogue of the sea anemone peptide ShK^'^ has 
recently entered clinical trial for autoimmune diseases'\ so the thera- 
peutic prospects for a potent, selective and stable peptide such as 
HsTXl [R14A] are quite promising. Analogues of other scorpion 
toxins also show promise in this respect, for example an engineered 
analogue of the three-disulfide-containing peptide BmKTX'^. 
The finding that R14 analogues of HsTXl preferentially inhibit 
the proliferation of CCR7" T^m lymphocytes and NA-NK lympho- 
cytes, leaving CCR7^ naive/TcM and A-NK cells able to perform 
their normal functions, suggests that our HsTXl analogues do 
not affect KCa3.1 channels and will not induce generalized 
immunosuppression. 

The novel HsTXl analogues described here display a remarkable 
selectivity for Kvl.3 over Kvl.l channels. Homotetramers of the 
Kvl.3 channels have a discrete tissue distribution, with expression 
in lymphocytes and the olfactory bulb'^. In other tissues expressing 
the Kvl.3 subunit, the resulting channel is a heterotetramer of 



different KvLx subunits. Selective blockers of Kvl.3 homotetramers 
are likely not to affect heterotetrameric channels and are expected to 
display no or few side effects, as was observed with the ShK analogue 
ShK-186^'^''\ In addition, as similar peptides do not cross the blood- 
brain barrier, HsTXl and its analogues are unlikely to reach and 
affect Kvl.3 homotetramers in this tissue. 

Our results also provide an excellent illustration of the predictive 
power of MD simulations and umbrella sampling PMF calculations, 
and how they can be exploited to design peptide analogues with 
enhanced selectivity properties''^'''* Because experimental structures 
for protein-ligand complexes are lacking in most cases, accurate 
construction of a complex structure using docking and MD simula- 
tions is of critical importance as it provides a basis for engineering 
modified peptides. In the next step, the free energy change associated 
with a given mutation can be accurately calculated using either the 
PMF method or free energy perturbation if the binding mode is 
preserved. In HsTXl, comparison of its binding modes with Kvl.l 
and Kvl.3 almost uniquely identified R14A as a mutation that should 
improve the Kvl.3/Kvl.l selectivity, and this was supported by the 
calculated binding free energies and confirmed by the measured 
binding constants. 
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Methods 

Modelling and docking. Here we provide a brief description of the computational 
methods and refer to our previous work^^-^'' for details. The structure of 
HsTXl[R14A] was generated from that of HsTXl using the mutator plugin in the 
VMD software^". Because the affinity of HsTXl for Kvl.2 is already below the 
measurement threshold, and R14A would reduce the toxin's affinity for Kvl .2 further, 
we consider here only the complexes with Kvl.l and Kvl. 3. 

For the Kvl.l and Kvl. 3 structures, we used homology models based on the Kvl.2 
structure, which were developed in^"^. The initial poses for the Kvl.x-HsTXl com- 
plexes were found using the docking program HADDOCK^\ and then refined in MD 
simulations. The most stable complex in each case was further run for 20 ns to check 
equilibration and for data analysis. MD simulations were performed using the 
NAMD program^' with the CHARMM36 force field^l The PMF of each toxin was 
determined from the umbrella sampling MD simulations using the weighted his- 
togram analysis method^^'^'*. 

The binding constant was determined from the integral of the PMF, W{z), along 
the reaction coordinate using equation 1 

Ke, = rf Te-^/'^^t/z (1) 

Jzl 

Here Zx was in the binding pocket and in the bulk where W vanishes. The factor nB? 
measures the cross sectional area of the binding pocket. The standard binding free 
energy of the toxin was obtained from equation 2 

G^=-/:nn(K,gCo) (2) 

where Co is the standard concentration of 1 M. 

Peptide synthesis. HsTXl was synthesized on a Prelude peptide synthesizer using an 
Fmoc-tBu strategy. The peptide was synthesized starting with Rink amide resin 
(Peptides International, Louisville, KY). Prior to coupling the Arg at position 14, the 
resin was split into equal portions replaced into separate reaction vessels and coupled 
with either Arg, Ala, Phe, Abu, or Val. The remainder of the sequence was 
subsequently continued until the full length peptide was assembled. All couplings 
were mediated with diisopropyl carbodiimide and 6-chloro-hydroxybenzotriazole. 
Following solid-phase assembly of the linear peptide chain, the peptide was cleaved 
from the solid support and simultaneously deprotected using Reagent K for 2 h at 
room temperature. The crude peptide was precipitated into ice cold diethyl ether and 
washed thoroughly to remove cationic scavengers from the cleavage cocktail, 
dissolved in 50% aqueous acetic acid, then diluted in water and the pH adjusted to 8.0 
with NH4OH. Disulfide bond formation was facilitated with reduced and oxidized 
glutathione according to previously used protocols for ShK^°'^^. 

The progress of folding was followed by RP-HPLC using a Phenomenex Luna CI 8 
column using a gradient of acetonitrile versus H2O containing 0.05% TFA from 10- 
70% over 35 min. Folding of the three disulfide bonds was also confirmed by the loss 
of 6 mass units from the crude material as determined by ESI-MS. 

NMR spectroscopy. Samples were prepared by dissolving lyophilised HsTXl and 
HsTXl[R14A] in 90% H2O/10%'H2O, pH 4.0, to concentrations of 500 and 400 |aM, 
respectively. One-dimensionaPH spectra and two dimensional homonuclear TOCSY 
spectra with a spin lock time of 70 ms were acquired at 35°C on a Bruker Avance 
600 MHz spectrometer. A NOESY spectrum (mixing times 200 ms) was also 
acquired at 35°C. All spectra were processed in TOPSPIN (version 3.2, Bruker 
Biospin) and analysed using CcpNmr- Analysis (version 2.1.5). ^H chemical shifts 
were referenced to the residual water signal. Chemical shift assignments for backbone 
and side chain protons of HsTXl [R14A] were made by conventional analysis of 
TOCSY and NOESY spectra. A complete assignment of the^H NMR resonances of 
HsTXl [R14A] was obtained. 

CD spectroscopy. HsTXl and HsTXl [R14A] were dissolved in 10 mM acetate buffer 
at pH 3.8 and CD spectra were recorded on a Chirascan v4.2. 12 spectropolarimeter at 
concentrations of 94 and 89 |aM, respectively, at 20°C. Spectra were collected at 1 nm 
intervals over the wavelength range 260-190 nm in a 1 mm path-length cuvette. 

Cells and cell lines. Mouse L929 fibroblasts stably transfected with mKvl.l or 
mKvl .3^^ were kind gifts from Dr. K. George Chandy (University of California, Irvine, 
USA). 

CD45RC"CCR7" ovalbumin- specific Ova-GFP Tem lymphocytes were derived 
from Lewis rats by Dr. Fliigel (University Medical Center Gottingen, Germany)^^. 
They were maintained in culture through alternating rounds of expansion in cyto- 
kine-rich medium and antigen-induced activation; as antigen-presenting cells for the 
Ova-GFP Tem lymphocytes, we used irradiated Lewis rat thymocytes'^. Lewis rat 
splenocytes were isolated using Histopaque-1077 gradients from the thymus 
donors'*^. All procedures involving vertebrates were conducted under a protocol 
approved by the Institutional Animal Care and use Committee at Baylor College of 
Medicine. Female Lewis rats (7-9 weeks old) were purchased from Charles River 
Laboratories (Hollister, CA) and housed with food and water ad libitum in a facility 
approved by the Association for Assessment and Accreditation of laboratory Animal 
Care International (AAALAC). Rats were euthanized by deep anesthesia foUowed by 
cardiac puncture. Death was ensured by decapitation. 



Buffy coats were purchased from the Gulf Coast Regional Blood Center (Houston, 
TX) and NK lymphocytes were further isolated by negative selection using the 
RosetteSep kit (Stemcell Technologies, Vancouver, Canada)'^. Cells were incubated 
overnight in medium supplemented with 500 lU/ml IL-2 and 1 ng/ml IL-15 for 
separation of adherent and non-adherent NK cells^^. The Institutional Review Board 
at Baylor College of Medicine determined that this work does not constitute human 
subject research as all buffy coats were de-identified and coded to preclude tracing a 
sample back to its donor. 

Electrophysiological analysis. Patch-clamp experiments were conducted in the 
whole-cell configuration using either a manual setup or a Nanion Port-a-Patch setup, 
both connected to an EPCIO-USB amplifier^'^l Chip or pipet resistances averaged 
2 MQ and the holding potential was set to — 80 mV. Series resistance compensation 
of 80% was used when currents exceeded 2 nA. Kv currents were elicited by repeated 
200 ms pulses from —80 to 40 mV, applied every 30 s, in normal Ringer solution 
with a calcium-free pipette buffer containing (in mM): 145 KF, 10 HEPES, 10 EGTA, 
2 MgCl2, pH7.2, 300 mOsm. IC50 values were determined by fitting the Hill equation 
to the reduction of area under the curve measured at 40 mV. 

Lymphocyte proliferation assays. Effects of HsTX-1, HsTX-l[R14A], and ShK-186 
on the proliferation of lymphocytes were measured by [^H] thymidine incorporation 
assay in 96-well plates, as previously described^*^'^°. Rat splenic mononuclear cells were 
plated at 10^ cells/well; Ova-GFP Tem cells at 50,000 cells/ well in the presence of 
irradiated thymocytes (2 X 10*^/well); and adherent or non-adherent natural killer 
cells at 50,000 cells/well. The Kvl. 3 blocking peptides were incubated with the cells 
for 30 min at 37°C before addition of the stimuli ( 1 |ag/mL concanavalin A for splenic 
T cells, 10 |ig/mL ovalbumin for Ova-GFP Tem cells, and 40 ng/mL PMA + 500 nM 
ionomycin for NK cells). Cells were incubated for a total of 72 h at 37°C and [^H] 
thymidine was added to the wells 16-18 h before harvesting of DNA on glass filters 
and counting of incorporated [^H] thymidine using a (3-scintillation counter 
(Beckman Coulter, Brea, CA). 

Proteolysis assays. Proteolysis assays were performed at a 250 : 1 substrate (peptide)/ 
enzyme ratio with pepsin, trypsin and a-chymotrypsin. For all assays, peptides were 
incubated with enzyme at 37°C for 4 h. As a positive control to ensure that active 
enzyme was present, bovine serum albumin was used as substrate in place of 
HsTXl [R14A]. All digestion assay data was analysed by reversed-phase HPLC 
(5-100% acetonitrile gradient, 10 min). Trypsin (EC 3.4.21.4, Sigma) and 
a-chymotrypsin (EC 3.4.21.1, Sigma), stocks were prepared in 50 mM Tris, 100 mM 
NaCl buffer (pH 7.4), and pepsin (EC 3.4.23.1, Sigma) stocks were prepared in 
10 mM acetic acid/HCl buffer (pH 2). The trypsin and a-chymotrypsin reactions 
were halted with 2.5% acetic acid (25% v/v), and the pepsin reactions were halted with 
2.5% of 200 mM glycine-NaOH buffer (pH 11). 

Haemolysis assays. Erythrocytes were isolated from buffy coats obtained from the 
Gulf Coast Regional Blood Center (Houston, TX), as described previously^\ Isolated 
haematocrit was diluted in PBS and incubated for 24 h at 37°C with 100 nM of the 
experimental peptides in 96-well plates. Absorbance of cell supernatants was then 
measured at 450 nm using a Powerwave XS microplate spectrophotometer (Biotek 
Instruments, Winooski, VT, USA). 
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